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We present a new method for computing the best approximation to a Killing vector on closed 
2-surfaces that are topologically 5*^. When solutions of Killing's equation do not exist, this method 
is shown to yield results superior to those produced by existing methods. In addition, this method 
appears to provide a new tool for studying the horizon geometry of distorted black holes. 
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INTRODUCTION 

An exact geometric sphere possesses a two parameter 
family of rotational Killing vectors while even a slightly 
distorted sphere may possess no Killing vectors what- 
soever. Nevertheless, one could imagine defining per- 
turbations of these initial vectors which would, in some 
well-defined "best" sense, represent the closest available 
approximation to vectors which almost satisfy Killing's 
equation on such a slightly distorted sphere. In this paper 
we introduce a definition which is best in a least squared 
sense and discuss some of its attributes. 

In general relativity, rotational Killing vectors play an 
important role in providing a quasi-local definition for 
the spin of a rotating body. A system of astrophysical 
interest, such as a pair of orbiting black holes, possesses 
no global rotational Killing vectors. In this case, angular 
momentum can only be rigorously defined for the sys- 
tem as a whole in terms of asymptotic rotational Killing 
vectors. 

A quantity of great importance in the evolution of 
black-hole binaries is the spin of the individual black 
holes. The spin of such black holes can only be deter- 
mined by some approximate quasi-local definition (see 
Ref. [H for a review). There exist many different quasi- 
local definitions for angular momentum, but they all take 
the form of an integral over a 2-surface with topology S'^, 
and all require a rotational Killing vector on this surface. 

In numerical relativity, the quasi-local spin of a black 
hole is most often expressed as 



'(C) 



(1) 



where Kij is the extrinsic curvature of a spatial hyper- 
surface with metric 7ij, d?S^ is the area element of an 
surface of integration S taken to be a black hole's 
apparent horizon, and is a Killing vector of the met- 
ric hij induced on S by ^ij. This form was derived by 
Brown and YorkQ and later within the Isolated Horizons 
framework (see Ref. 0] for a review) . It gives the angular 
momentum of the rotation associated with the rotational 
Killing vector Unfortimately, the induced metric hij 



will not admit a solution of Killing's equation for the case 
of orbiting black hole binaries. In this situation, one has 
no recourse but to find some reasonable approximation 
for the Killing vector required in Eq. ((T|). 

In some cases, conformal Killing vectors have been 
used, and have yielded physically reasonable results 4 1. 
Better still, a "Killing Transport" (KT) tcchnique[5|, 
which finds exact Killing vectors when they are present, 
has been recently adopted, and appears to give physi- 
cally reasonable results (cf Refs. 0, [a|). Our definition is 
shown to be even better in a well-defined least squared 
sense and, for coalescing binary black holes, yields other 
interesting results worthy of further investigation. 



EQUATIONS FOR THE BEST APPROXIMATE 
KILLING VECTOR 



An arbitrary vector field ^' on 5*^ can be decomposed 
into two scalars d and v 



(2) 



where Di is the covariant derivative compatible with the 
metric hij induced on the surface S and eij is the Levi- 
Civita tensor. Similarly, the general gradient of a 1-form 
can be expressed as 



DiS^j = Lei 



h,jA 



(3) 



where L and A are scalars, and Sij is symmetric and 
trace-free. Eqs. ([2]) and ^ immediately imply that 



A 



L = -\D'D,v, 
5y = D,D,d-\h,jD''Dkd 



(4) 
(5) 
(6) 



For ^* to be Killing, it must satisfy Killing's equation 
^(i^j) ~ ^ where the parentheses denote symmetrization. 
This implies that A and Sij must vanish if ^* is Killing. 
We may choose A to vanish, in which case Eq. (j4|) implies 
that d is harmonic. But, assuming a non-singular metric, 



2 



the only harmonic function on S'^ is a constant, which 
makes no contribution to Thus, we arc left with 

C ^ ^'D,v, (7) 
ACj = Le,j+S,j. (8) 

Our goal is to find an approximate Killing vector that 
minimizes the non-Killing aspects of Clearly, having 
already set A to zero, a solution with Sij as close to zero 
as possible is what we desire, so we proceed by finding a 
vector that minimizes SijS^^ . 

From Eqs. ^ and © it follows that 

L^^e,,D'e, (9) 

and 

S,,S'' = {D,Djv){D'D^v) - ^iD'^Dkvf. (10) 

Now, we wish to cxtremize SijS^^ with respect to v. But, 
we must do this in a way that is independent of the nor- 
malization of Using l^p = {Div){D'^v), we choose the 
following scalar function on S 

C = S,,S'' + \^m{Dkv){D''v), (11) 

where is the Ricci scalar associated with hij and 
is a dimensionless constant. Varying with respect to w, 
SC/Sv = yields a fourth-order scalar elliptic equation 
for V that can be rewritten as a pair of second-order scalar 
elliptic equations for v and L 

D'DiL-{l-e)[^iD'^R)D,v-^RL] =0, (12) 

DW,v + 2L = 0. (13) 

Eqs. ([T^ and (fT5|) can be solved for L and v, with 
the Lagrange multiplier O fixed by the requirement that 
Eq. be intcgrable on S'^. Given a solution for v, the 
approximate rotational Killing vector is given by Eq. ([7]) 
once it is normalized to have affine length 27r. We note 
that Eq. p2|) is satisfied by a Killing vector ^* when = 
(with L defined by Eq^ © and D,v = i2/'^R)DiL). 

It has been noted [7|, |8| that an approximate Killing 
vector should be divergence free D^^* = 0, in part be- 
cause the angular momenta computed using such an ap- 
proximate Killing vector possess a certain gauge invari- 
ance. Approximate Killing vectors constructed using the 
KT method are not guaranteed to be divergenceless, al- 
though at least in certain cases Q this can be enforced a 
posteriori. These approximate Killing vectors also inherit 
an additional problem. A defining equation of the KT 
method which can be written in the form of Eq. ([9]), 
is not satisfied by solutions of the KT equations unless 
there is a Killing vector. This is due to the path depen- 
dence of the solution scheme. Moreover, Eq. ^ cannot 
be enforced a posteriori. 

Because our approximate Killing vector is defined from 
the solution of Eqs. (fT2|) and (fTS]) via Eq. ([7]), it is guar- 
anteed to be divergenceless. Furthermore, because our 
solution is obtained via a global solution of elliptic equa- 
tions, Eq. ^ is also guaranteed to be satisfied. 



TESTS 

We have implemented a code to solve Eqs. and 
(|13p for situations where the metric hij is conformal to a 
unit 2-sphere: 

ds^ = V'V2(d6|2 + sin^ Odcp^). (14) 

Details of the solution scheme will be presented in a fu- 
ture paper 0. While this form for the metric may seem to 
be a strong simplification, a scheme based on this form is 
suitably general since any sufficiently smooth metric on 
5^ is conformally equivalent to a unit 2-sphcrc. 

To explore our method, we first consider the case where 
the conformal factor is given by an £ = 2. to = scalar 
spherical harmonic with its axis of symmetry rotated to 
a direction given by {9', (j)') 

2 

i;{9, c^)=A^B Y2m{9, (f)Y*^{9', 0'), (15) 

m=-2 

and A and B are real constants chosen to guarantee that 
V' > everywhere. The form of Eq. (fT5|) guarantees that 
the metric possesses a rotational Killing vector and in 
all the cases we attempted, solving Eqs. ([T^ and (fT^ 
returned a solution where the axis of symmetry was cor- 
rectly rotated by {9' , cj)') and for which = 0. As men- 
tioned above, the solutions are divergenceless and satisfy 
Eq. Q to the level of roundoff error. Furthermore, we 
find SijS"^^ = also to the level of roundoff error, as 
expected for a solution that yields a true Killing vector. 

Interestingly, Eqs. (|12p and (|13p in general have multi- 
ple solutions. In fact for the example given by Eq. (jlSp . 
there exist an infinitely degenerate set of solutions where 
the axis of approximate symmetry lies anywhere in the 
rotated equatorial plane. For these solutions, 7^ and 
SijS^^ ^ since these solutions are not true Killing vec- 
tors. Again, the solutions are divergenceless and satisfy 
Eq. ([9]) to the level of roundoff error. 

We have also tested the system of equations against 
the numerically generated initial data for corotating and 
non-spinning equal-mass black-hole binaries as described 
in Ref.[l]. We use data for different orbital separations, 
parameterized by the dimensionless orbital angular ve- 
locity MHq, where M is the total irreducible mass of 
the binary. For the case of corotating black holes, the 
spin of each black hole is aligned with the direction of 
the orbital angular momentum. In all cases tested, we 
have found that the measured spins based on an approx- 
imate Killing vector obtained using Eqs. and are 
nearly identical to those based on the KT method, with 
differences growing only to a few parts in 10^. While 
the resulting spins are nearly identical, our solutions are 
measurably different from the results of the KT method, 
and we find that our system of equations produces so- 
lutions for which ) = (47r)~^ ^ SijS'^^dil is always 
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FIG. 1: The value of {SijS''-') for the new method is dis- 
played as a sohd (black) line and its value when using the KT 
method is shown as a dashed(red) line. The analogous quan- 
tity (2A^) when using conformal Killing vectors is shown as 
a dot-dot-dashed(black) line visible in the inset where a loga- 
rithmic scaling is used. The value of the Lagrange multiplier 
Q obtained using the new method is also displayed as a dot- 
dashed(blue) line and its scale is shown on the right. All data 
is from one corotating black hole in an equal-mass binary. 
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FIG. 2: The top portion of the figure displays the directions 
of the additional approximate symmetry axes with </!> = in 
the positive x direction. The axis direction of the Solution 1 
is displayed as a solid(black) curve and the direction for So- 
lution 2, with tt/2 subtracted, is displayed as a dashed(red) 
line. The lines are visually coincident. The lower portion of 
the figure displays At/i = 01 — 02+^/2 showing the degree to 
which the two axes are approximately orthogonal. All data is 
from one corotating black hole in an equal-mass binary. 



smaller than that produced by the KT method. Note 
that in both cases, the approximate Killing vectors have 
been properly normalized to have afline length 2tt. Fig- 
ure[T]shows (SijS^^) as computed by our method and the 
KT method, and the analogous quantity (2A^) obtained 
by using conformal Killing vectors. On a separate scale. 
Fig. [T] also shows the value of O obtained by our method 
for the same data. 

As with the first example given by Eq. (jlSp . our equa- 



tions yield multiple solutions for the numerically gen- 
erated black-hole data. The solutions shown in Fig. [T] 
have their axis of approximate symmetry aligned with 
the orbital angular momentum (z-axis) as expected. In 
all cases, we find two additional solutions. Both have 
their axis of approximate symmetry in the orbital plane 
and we refer to them as Solutions 1 and 2. We find 
that for all orbital separations, the directions of these 
two symmetry axes are orthogonal to the level of nu- 
merical error in the code (truncation error). For large 
separation (small Mflo), Solution 1 has its axis of ap- 
proximate symmetry pointed roughly in the direction of 
motion of the black hole (y-axis) at an azimuthal angle 
of , and Solution 2 has its axis pointed roughly toward 
the companion black hole (x-axis) at an angle of 02 ■ The 
directions 4>i and (f)2 (mod tt) for these two solutions are 
displayed in Fig. [2] Here, zero azimuthal angle is in the 
direction of the positive x-axis. For small separations, 
the roles are reversed and we find that Solution 1 points 
roughly toward the x-axis and Solution 2 toward the y 
axis. 

The regime between large and small orbital separa- 
tions, where the two solutions swap orientations, is quite 
interesting. Over most of the range of Mfio considered, 
4>i and 4'2 change gradually. However, over a narrow 
range of M^lg, the angles change rapidly but smoothly, 
rotating by an angle of approximately tt/2. Interestingly, 
the value of {SijS'^^) for Solution 1 is smaller than that 
for Solution 2 for all separations. At the point where 
4>i = 02 — 7r/2 = 7r/4, curves of {SijS'^^) for the two solu- 
tions "appear" to cross. However, a careful examination 
shows this to be an "avoided" crossing. We shall examine 
this behavior in more detail in a future paperQ. 

Finally, if we measure the spin of the black holes us- 
ing the approximate Killing vectors associate with So- 
lutions 1 and 2, all cases yield zero to truncation error. 
So, for the case of corotating equal-mass black hole bina- 
ries, we find three "orthogonal" solutions and only one 
of them yields a non- vanishing spin. 

We find similar results for the case of non-spinning 
black-hole-binary initial data. There is a solution with 
an axis of approximate symmetry aligned with the or- 
bital angular momentum and two additional solutions 
with their axes in the orbital plane and similarly orthog- 
onal to each other. As discussed in Ref. [1], the non- 
spinning black hole data we use are defined by setting 
the spin measured by the KT method to zero. The cor- 
responding spins measured by our method again differ at 
most by a few parts in 10'', and in all cases we find the 
value of {SijS'^^) for the new method to be smaller than 
the value from the KT solution. The behavior of the 
additional solutions is qualitatively the same as seen for 
the case of corotation and is displayed in Fig. [3l While 
the behaviors are generally similar, the directions of the 
approximate symmetry axes are somewhat different and 
the rapid change in the direction of the solutions occurs 
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FIG. 3: The directions of the additional approximate symme- 
try axes from one non-spinning black hole in an equal-mass 
binary. See description for Fig. [2] 



at much smaller separation. 



DISCUSSION 

We have mentioned two previous methods defined in 
the literature for use in computing the spin of rotat- 
ing black holes that lack axial symmetry. Both have 
been considered useful in the past, and yet both have 
shortcomings. When a Killing vector does not exist, the 
conformal Killing approach returns a vector for which 
DiC 7^ 0- By contrast, the Killing Transport method 
constructs a ^' that can often be made divergenceless. It 
also constructs the scalar L (see Eq. ([3])), but generally 
this does not satisfy Eq. Our new method not only 
ensures both that is divergenceless and that Eq. ([9]) 
is satisfied, but it also is best in the sense that (SijS^^) 
is minimal. Since a Killing vector cannot be produced 
where one does not exist, the usefulness of our results for 
an approximate Killing vector will depend on the extent 
to which physical questions (such as concern black hole 
spins) can be given meaningful answers. In particular, it 
may have immediate application in giving a more refined 
definition for binaries containing black holes without in- 
dividual spin. 

One sense in which our results already appear mean- 
ingful relates to the vectors found to reside in our xy- 
plane, for which the corresponding spins are computed to 
be zero to the level of truncation error. Any other out- 
come for these would have been somewhat unpalatable. 
Their actual orientation for binary black holes at large 
separation can be easily interpreted in terms of boosted 
frames. Closer in, their combined dramatic rotation at 
some critical separation warrants further investigation, 
as does the apparent occurrence of avoided crossings for 
{SijS'^^) and for O at the critical separation. 

Another sense in which our results appear meaningful 



relates to the solution associated with our z-direction. 
For sufficiently large black hole separation, our result 
leads to spins which are in good agreement with those of 
the KT method. This is not unreasonable, even though 
we find small differences between the two solutions for 
^* over the surface of the apparent horizon. However, 
for highly distorted black holes — such as near merger 
or with near maximal rotation, or for unaligned spins — 
we can imagine that our results could prove to be more 
robust. More extensive comparisons than we have been 
able to carry out here will be necessary before this ex- 
pectation might be practically substantiated. 

To a high degree of precision, it appears that the axes 
of the approximate symmetries we find form an orthog- 
onal basis in all the cases we have examined. This is 
of considerable interest, because this basis is directly re- 
lated to intrinsic properties of the apparent horizons we 
have studied. Thus, in addition to the reasonableness of 
our results associated with both the rotation axis and the 
orbital plane, our method appears to give a new tool for 
studying the horizon geometry of distorted black holes 
in general, and the volatile dynamics of black holes dur- 
ing collision in particular. Further investigation of the 
usefulness of this new tool will be forthcoming. 
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